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We propose a bigravity analogue of the F(R) gravity. Our construction is based on recent ghost- 
" free massive bigravity where additional scalar fields are added and the corresponding conformal 

transformation is implemented. It turns out that F(R) bigravity is easier to formulate in terms 
. of the auxiliary scalars as the explicit presentation in terms of F(R) is quite cumbersome. The 

£S| 1 consistent cosmological reconstruction scheme of F(R) bigravity is developed in detail, showing the 

possibility to realize nearly arbitrary physical universe evolution with consistent solution for second 
metric. The examples of accelerating universe which includes phantom, quintessence and ACDM 
acceleration are worked out in detail and their physical properties are briefly discussed. 
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I. INTRODUCTION 



The formulation of massive spin-two field or massive graviton has a long history initiated from the free field 
_ formulation by Fierz and Pauli |l| (for recent review, see ]2|). In spite of the success of the free theory, it has 
. been known that there appears the Boulware-Dcser ghost [3[ in the naive non-linear extension of the Fierz-Pauli 
formulation. Furthermore, it has been also known that there appears a discontinuity in the limit of m — > in the 
free massive gravity compared with the Einstein gravity. This discontinuity is due to the extra degrees of freedom 
in the limit and is called vDVZ (van Dam, Veltman, and Zakharov) discontinuity Q. The extra degrees of freedom 
^0 \ can be screened by the non-linearity, which becomes strong when m is small. Such mechanism is called the Vainstein 
mechanism A similar mechanism works [(| for the bending mode of the so-called DGP model Q- Moreover, the 
' scalar field models, where the Vainshtein mechanism works, have been proposed. 

Recently, there has been much progress in the non-linear formulation of the massive gravity 0, Q without the 
f**«« 1 Boulware-Deser ghost Although the corresponding formulation of massive spin- two field is given in the fixed or 
non-dynamical background metric, the ghost-free model with the dynamical metric has been also proposed (for 
the recent cosmological aspects of massive ghost-free and bigravity models, see Since the corresponding model 

contains two kinds of symmetric tensor fields, the model is called bi- metric gravity or bigravity. The massive gravity 
^ ■ was applied in Ref. |12{ to explain the current accelerating expansion of the universe. The accelerating cosmology in 
terms of the recent formulation of the ghost-free bigravity was discussed in [l3[ . 

It is commonly accepted nowdays that the expansion of the current universe is accelerating. This was confirmed 
by the observation of the type la supernovae at the end of the last century In order that the current cosmic 

acceleration could occur in the Einstein gravity, we need the mysterious cosmological fluid with the negative pressure 
called dark energy (for recent review, see [la|). The simplest ACDM model of dark energy is composed of the 
cosmological term and CDM (cold dark matter) in the Einstein gravity. The ACDM model, however, suffers from the 
so-called fine-tuning problem and/or coincidence problem. In order to avoid these problems, many kinds of dynamical 
models have been proposed. 

Among such dynamical models, much attention has been given to the so-called F(R) gravity which was proposed 
as gravitational alternative for cosmic acceleration in Refs. [It], HH (for recent review, see [IB]). In F(R) gravity, the 
scalar curvature R in the Einstcin-Hilbcrt action is replaced by an appropriate function F(R) of the scalar curvature. 
In this paper, we propose a bigravity analogue of the F(R) gravity. We formulate the theory which respects the 
desirable properties of the recent bigravity models and for example, the Boulware-Deser ghost does not appear. It is 
demonstrated that the obtained field equations are consistent with each other and consistent cosmological solutions 
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can be obtained. Furthermore, we show that a wide class of the cosmological solutions, including the accelerated 
expanding universe, can be realized in this formulation. Therefore, the models under consideration have much richer 
structure than simple bigravity recently investigated in [l3| . 



II. GHOST-FREE F(R) BIGRAVITY 

A model of bi- metric gravity, which includes two metric tensors g^ v and was proposed in Ref. [Toj . The model 
describes the masslcss spin-two field, corresponding to graviton, and massive spin-two field. It has been shown that 
the Boulware-Deser ghost Q does not appear in such a theory. 

The action is given by 



S hi = M 2 g J d 4 x ^ ~ dot gR,^ + M 2 J d 4 x y / - det / R.W 

+2m 2 M 2 ff / d 4 xy/-detgJ2Pn e n (V^F 1 /) ■ (1) 

•J ^— n 



Here R^ is the scalar curvature for g^ and R^' is the scalar curvature for f^ v . The tensor \J g 1 / is defined by 
the square root of g pp f pv , that is, (y/ g" 1 /^) (\/g~ x f) = 9 tlp 'fpv For the tensor X^, e„(A)'s are defined by 



e (A) = l, ei (X) = [X], e 2 (A) = i([A] 2 -[A 2 ]), 
e 3 (A)-i([A] 3 -3[A][A 2 ]+2[A 3 ]), 

e 4 (A) = ^([A] 4 - 6[A] 2 [A 2 ] + 3[A 2 ] 2 + 8[A][A 3 ] - 6[A 4 ]) , 

e k {X) = for k > 4. (2) 

Here [X] expresses the trace of X: [X] = X 11 

We now construct a bigravity model analogous to the F(R) gravity. Before going to the explicit construction, one 
may review the scalar-tensor description of the usual F(R) gravity [la ]. In F(R) gravity, the scalar curvature R in 
the Einstcin-Hilbert action 

SeH = [ d 4 Xy/^g ( — + £ ma ttcr ) , (3) 



k 2k 2 

is replaced by an appropriate function of the scalar curvature: 

// p(R\ \ 
d 4 Xy/^g ( -y-^- + ^matter J ■ (4) 

One can also rewrite F(R) gravity in the scalar-tensor form. By introducing the auxiliary field A, the action Q of 
the F(R) gravity is rewritten in the following form: 

S = h- I ^^{F'iA) + F(A)} . (5) 

By the variation of A, one obtains A — R. Substituting A = R into the action ([5]), one can reproduce the action in 
((4]). Furthermore, we rescale the metric in the following way (conformal transformation): 

g^^e a g^, <j = -\nF'(A). (6) 

Thus, the Einstein frame action is obtained: 

V(a) = c°g (e-) - e 2 */ (g (e-)) = ^ - ^ . (7) 

Here g(e~ a ) is given by solving the equation a = — In (1 + f'(A)) = — lnF'(A) as A = g(c~ <T ). Due to the scale 
transformation the scalar field a couples usual matter. 
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In order to construct a model analogous to the F(R) gravity, we added the following action to the action ([T]): 

Si = -M 2 J dW-detg V^g>* u d^d v <p + V(ip)\ + J d i xC mat t cl - (e%„, $,•) . (8) 

Here we denote the matter field by <£v As discussed in [13], the action (j8]) does not break the good properties like 
the absence of the Boulwarc-Dcser ghost. 

By the conformal transformation — > eT^g^, the total action /Stotai = Shi + S\ is transformed to 

Stotal — > 

See. = M 2 f [ d 4 x^drtf R^ + 2m 2 M 2 s / d^v^det^^ Z?^*" 2 )^ (vF 1 ?) 

+M 2 j d^^A^g {e-fRM + e~ 2ifi V(^ + j d 4 xC mattei (g^, $0 • (9) 



Then the kinetic term of ip and the coupling of ip with matter disappear. By the variation over tp, we obtain 

4 

= 2m 2 M 2 ff ]T/3„ (| - 2) c(t- 2 )^e„ (VF 1 ?) + M fl 2 {-e"^ - 2e- 2 *VM + e^V (if)} . (10) 



n=0 



Eq. (fTU)) can be solved algebraically with respect to ip as if = ip^R^ 9 \e n y^J g^ 1 J^jJ . Then by substituting the 
expression of ip into a model analogous to the F(R) gravity follows: 

S FR = M) J d 4 aV"det/ R u} 

+M 2 J d 4 x^/-detgF e„ (V^F 1 /)) + / d*x£ ma . tteI (g^, *<) , 



Note that it is difficult to solve (|TU|) with respect to ^ explicitly. Therefore, it might be better to define the model by 

introducing the auxiliary scalar field ip as in Of course, in some cases can be explicitly 

found. For instance, in the minimal case, where (3q = 3, /?i = — 1, P2 = /?3 — 0, and /34 = 1, one may consider the 
simplest case V = Vc~ v with a constant Vq. Then Eq. (fTU)) reduces to 



= m 2 M 2 ff f-12e" 2v ' 



5 e (VF 1 ?) + 3e"^ei {VFl)) - M 2 e~* (flW + , 
which can be solved with respect to e _ 5 as 



(12) 



e 2 



eo ( VFV) ^ ^ 64e ( A /F T ?) * 12m 2 M 2 ff e^ eo ( ^=1/) 



(13) 



and we obtain 



F{R^\e n [^f) 



( 



^(VF 1 /) 



± 



(vfv)' 



M^e-v (R(s) + V ) 



8e ( v /<r T 7) ^ 64e (VF 1 /) 2 12m 2 M 2 ff e-^e (v^FV) 
/ 



[V~9 1 f) M$e-? (Mb) + V ) 



61 (^) ± 
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xV 



( ( 

-2 In 



V 



± 



V 



ex (V.9- 1 /) 






+ ^o) 


64e [Vg-'f) 


12m 2 M c 2 ff e- 


- 2 ^e 





(14) 



Hence, we may define the analogue of the F(R) gravity by ©■ 

Even for the sector including f^ v , one may consider the analogue of the F(R) gravity by adding the action of 
another scalar field £ as follows: 



By the conformal transformation for / M „: / M , y — > c - ^ f^ v , instead of ([9]), wc obtain 
S F = Mj J cPxy/-dctf {e^i?^+e" 2S C/(0} 

+2m 2 Af c 2 ff / d 4 x^d^J2P^~ 2)v ~^en(^/g~ T f) 

n=0 

+M 2 J d 4 x^-detg {e~ v R^ + e~^V{>p)] + J d 4 z£ ma ttcr (g^, $0 
Again the kinetic term of £ vanishes and by the variation of tp and £, we obtain 



(15) 



(16) 



= 



4 

2m 2 Af 2 ff ^/3 Il (^-2)e(^- 2 )^-t« e „( v ^ T 7) + M 2 {-e^i^ - 2e~^V&) + e" 2 ^'^)} , (17) 

n=0 

= -2m 2 Af 2 ff J2 ^e(i- 2 )^t? e „ (VF 1 /) + M) {-e"^ - 2 C - 2 «[/(£) + c- 2 «£/'(0} . (18) 

The obtained equations (fT7|) and (fT8|) can be solved algebraically with respect to <p and £ as 93 = 
95 (i?( ff ),E^),e„ (V.^ 1 /)) and £ = £ (i? (9) , e„ fVff" 1 /))- Substituting the expression of v? and £ into (JTBJ, 
we obtain a model analogous to the F(i?) gravity: 



57 



Af 2 / d 4 x^dct7F^) (RW,R(f) ien (V^ 1 /)) 



-2m 2 Mi 



d 4 :. 



F 



71=0 

+M 2 / dSv/^det^F^) (i?(9),i?(/),e n (v / F T 7)) + / d A x£ matteI (g^, $<) 
.^.^(VF 1 ?))^) 



+c 



(19) 



We should again note that it is difficult to explicitly solve (flTj) and (fl~8| with respect to <p and £ and it might be 
better to define the model by introducing the auxiliary scalar fields tp and £ as in (1161) . 
Hence, we succeeded to obtain the bigravity analogue of the F(R) gravity. 



III. COSMOLOGICAL RECONSTRUCTION 

Wc now consider the minimal case, where 



S hi = M 2 J d 4 x^/-dctgR^ +M 2 f J d A x yj - det / R {f) 

+2m 2 M 2 s J d 4 x^-detg (s - tr V^F 1 / + dot v 7 ^ 1 /) . (20) 



In order to evaluate Sy/g 1 /, we consider two matrices M and N, which satisfy the relation M 2 = N. Since 
SMM + MSM = SN, we find 

SM = 5NM- 1 - MSMM' 1 . (21) 

By using (|2"Tj) iteratively, one obtains 

SM = SNM- 1 - MSMM- 1 = 8NM- 1 - MSNM~ 2 + M 2 SMM~ 2 = ^{-l) n M n SNM- n - 1 . (22) 

n=0 

Then by carefully considering the trace of Eq. (|2"T1) , we find 

tr SM = -tr (M^SN) . (23) 



For a while, we work in the Einstein frame action (|20[) with (JSJ) and (fTSj) but the contribution from the matter is 
neglected. Then by the variation of g pv , one obtains 



= Mg 



(^9^R (9) ~ R^j + m 2 M cS (3 - tr y/f*f) + f w (VFl) 



+ \ {^9 p °d p yd ai p + V(<p)j 9^ - . (24) 

On the variation of f pvi we obtain 

= Mj {^f^R U) ~ R^j + m 2 M ctf j/^ (3 - tr y/f*f) - (v 7 ^ 1 /) 3 pT /™ 

+ ^ (^f pa d P (,d a ti + U(0) V - \d»Zd v Z ■ (25) 

We now assume the FRW universes for the metrics g pv and f pv : 

3 33 3 

ds 2 g = 9^dx^dx v = -dt 2 + a(t) 2 ^ (da; 4 ) 2 , ds) = ^ f pv dx^dx u = -dt 2 + b(t) 2 ^ (dx 1 ) 2 . (26) 

Then the (t, t) component of gives 

= -3M 2 H 2 3m 2 M 2 s (l - ^ - ^ 2 - 1^) , (27) 

and components give 

b\ 3 . 9 1, 



= M 2 (2H + 3H 2 ) + 2m 2 M 2 s ( 1 - - J - -<j? + -V{<p) . (28) 



Here H = a/a. On the other hand, the (t,t) component of (|2~5j) gives 



= -3M^ 2 - m 2 M 2 s ( 1 ^) |C 2 - \U{ 
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and (i,j) components give 

= M) (2K + 3K 2 ) + 2m 2 M 2 ff (l - ^ - -J? + . (30) 



Here K = b/b. Hence, 



= 2M 2 H - m 2 M 2 s ( 1 - -) - j^ 2 , 



= -M 2 (2H + 6iJ 2 ) - 5m 2 M 2 s fl-^j- V(cp) , (32) 

= 2M}K + m 2 M 2 s (l - - ^ 2 , (33) 

= -M) (2K + 6A" 2 ) - m 2 M 2 s (3 - ^ - [/(£) . (34) 

One now redefines scalar fields as ip = ip(r]) and £ = £(£) and identify 77 and £ with the cosmological time t. Then we 
find 

^ = 2M 2 H - m 2 M 2 ff f 1 - b - V (35) 



Here 



V(t) = -M 2 (2H + 6iJ 2 ) - 5m 2 A/ 2 ff ^1 - ^ , (36) 

^ = 2M 2 K + m 2 M 2 s (l - ~) , (37) 

= -M) (2if + 6if 2 ) - m 2 M 2 s ^3 - ^ . (38) 

0,(77) = 3^'W 2 , V( V ) = V(<p(r 1 )) , ^(0 = 3r(C) 2 , U(0 = U(£(0)- (39) 

Then for arbitrary a(t) and 5(i), if we choose w(i), t^(i), c(i), and U(t) to satisfy Eqs. (|35H38p . a model admitting 
the given a(t) and 6(i) evolution can be reconstructed. 

Consider the possibility not to introduce the extra scalar field x fl5|) - Instead of the introduction of we assume 
the metric in the following form: 

3 3 

ds 2 = f^dx^dx" = -c{tfdt 2 + b(t) 2 Y (dx*) 2 ■ (40) 



Then instead of Eqs. (|27H30|) . one gets 







3Af 2 ff 2 - 3™ 2 Af 2 ff ( 1 - ±) - > - \ V { V ) , (41) 



= M 2 (2H + 3H 2 ) + m 2 M 2 s ( 3 - c -2-)-h 2 + \v{<p) , (42) 



2 

a) A r 2 

= -MI 2 f K 2 - m 2 M 2 s ^3 - 2c - ^ c 2 , (43) 
= M) (2K + 3K 2 - 2LK^j + m 2 M 2 B ^3 - c - — ^ c 2 . (44) 

Here L = c/c. 

For a given a = a(t), Eqs. (|43p and (|44|) could be solved with respect to & and c. On the other hand, as in (|35|) and 
, Eqs. (|4"Tj) and (|4"2")l can be rewritten as 

^ = 2M 2 H-m 2 M 2 Jc-^\ , (45) 
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V(t) = -M 2 g (zH + fiff 2 ) - m 2 M 2 s (Q-c-5 b -) ■ (46) 



Here u>(t) and V(t) are denned by (|30[) . Then for arbitrary a(t), if we choose u(t) and V(t) to satisfy Eqs. (1431) and 
(|4T>1) . a model admitting the given a(i) can be reconstructed. 



IV. EXAMPLES OF ACCELERATING COSMOLOGICAL SOLUTIONS 

Let us consider several examples. As discussed around (O, the physical metric, where the scalar field docs not 
directly coupled with matter, is given by multiplying the scalar held to the metric in the Einstein frame in ([5J or ([20]) : 



9%T = e%„ • (47) 
The metric of the FRW universe with flat spatial part is conformaly flat and therefore given by 

ds 2 = a{tf (^-dt 2 + (<*=*) J • (48) 

In case a{t) 2 = j?, the metric (|48|) corresponds to the de Sitter universe. On the other hand if d(t) 2 = ^ with 1, 
by redefining the time coordinate by 

l n 

dt = ±—dt, (49) 



that is, 



the metric (148 II can be rewritten as 



m 

t = ± -t 1 "", (50) 

n — 1 



ds 2 = -dt 2 + (±(n 1 " ( dxl ) 2 ■ ( 51 ) 

Then if < n < 1, the metric corresponds to the phantom universe, if n > 1 to the quintessence universe, and if 
?i < to decelerating universe (For similar scenario in usual non- linear massive gravity, see also (l9T|). In case of the 
phantom universe (0 < n < 1), one should choose + sign in ± of (|49|) or (|50)) and shift t as t — > < — to- Then t = t n 
corresponds to the Big Rip and the present time is 7; < to and 7j — > 00 to the infinite past (t — > —00). In case of the 
quintessence universe (n > 1), we may again choose + sign in ± of (|4^1) or ([5U1) . Then t — > corresponds to t — >• +00 
and t — > +00 to t — > 0, which may correspond to the Big Bang. In case of the decelerating universe (n < 0), we may 
choose — sign in ± of (|49l) or ([50)) . Then i — >• corresponds to t —> +00 and t —> +00 to t —> 0, which may correspond 
to the Big Bang, again. We should also note that in case of de Sitter universe (n = 1), t — > corresponds to ? — > +00 
and t — > ±00 to t — > —00. Let us now choose the metric in the Einstein frame to be flat, where H = 0, and 



l 2n 



(52) 



Using l[39|) . we find 



and Eq. (|35|) gives 



12n 2 . , 

w(t) = -p- , (53) 



6n 2 

&~1= 2Bf2 . 2 - (54) 



Eq. (|54|) shows the behavior of the metric / M „: 



fin 2 \ 2 

= ^" = ( 1 + ^M^J (55) 
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Then for large t, we find fij —> Sij, that is, the flat metric. On the other hand, for small t 

36n 4 



'eff 1 



which becomes larger and larger. Since small t corresponds to large physical time t for the phantom, the dc Sitter, 
and the quintessence universes, the late-time acceleration could be generated by the evolution of f^v. 
Using (J35J), the potential is 



We also find 



t- 



12nf 36n_ / i , 2n' 

7^1703 • m 2 M 2 t 4 l m 2 M 2 „t 2 



= „. ^ - = off - v ° / . (5 8 ) 

1 + 



With the help of d3T|) and <J38J), we obtain 

144M?n 2 



12n 



*(*) = f V TT 7 (59) 



1 

72ra 2 M 



6n 2 
in 2 M 2 ff t 2 



f 2 



! A/ 2 . f i 4 I, 1 + 



ud) = - - t- v 6 „; / + ^ ^ + ] . (60) 



,2>l..r2 72 



When t is small, a(t) behaves as 

/ 8M 2 f A l 

"M~(snfc- u "'Jir ( 61 ) 

In order to avoid the ghost, we require <r(t) > 0, which gives a constraint for the parameters as follows: 

2M? 

> in 2 . (62) 



m 2 M 2 ff 



On the other hand, when t is large, the second term dominates in Eq. ([5 



<K*) ^- • (63) 

Therefore, o"(i) becomes negative although there does not appear the Boulware-Deser ghost there could appear an 
additional ghost associated with the scalar field £. We should also note the negative a conflicts with (|39[) and therefore 
the model cannot be identified with the analogue of the F(R) gravity. This problem can be, however, avoided by 
modifying the large t behavior. Indeed, large t does not always mean the late-time when we choose the physical time 
t in ([50|) as discussed after Eq. (fSTj) . In case of the phantom universe (0<n<l),£— s-oo corresponds to the infinite 
past (t — > -co). In case of the quintessence universe (n > 1) or the decelerating universe (n < 0), the limit of t — > +oo 
corresponds to that of t — > 0. Even in case of de Sitter universe (n = 1), t — > ±oo corresponds to t — > — oo. Therefore, 
the modification of large t does not affect the late-time behavior of the universe. 
Finally, the ACDM-like universe may be reconstructed: 

ds 2 = -dP + A 2 sinh 3 l - V (dx l ) 2 . (64) 



I ^ 

i=l 



Here A and I are constants. Changing the time variable t by 



dt = M , (65) 
A sinh 2 1 



i 
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I gives 



we obtain the conformal form of the metric as in (|48| . Eq. ([65]) 

lV2 ( 21 3 1 N 
A \ '4' 2, 



Here B(a;, a, b) is the incomplete beta function defined by 

B(x, a, b) 

Turn 



dxx - 1 (1 - .t) 6 " 1 



3 _P 

'4' 2, 



c^ = a(t) 2 = A 2 sinh 3 , t = -1 In (iT 
Here -B -1 (y, a, 6) is the inverse function of B(x, a, b) defined by x = -B~ 4 (y, a, 6) for y = B(x, a, b). Eq. (j3"9")l gi 



(66) 



(67) 



(68) 



Therefore Eqs. (J35J and give 



sinh -V- cosh 2 . 
t 2 i I 



27 A2 sinh ±® cosh 2 t® 
2m 2 M e 2 ff /2 Smfi Z C ° sn / 

have used 



27A 2 

6 = 1 J 

Here we 



135 .4 2 / 



(*) T - r , \ 135 A 2 .,*(*) , 2 * (*) 
■7^- , 7(t) = -7777- sinh -y- cosh 2 -V 2 

l|6"5"]) and (f6"8"]) . Hence, we find 

^ sinhl g cosh g (cosh 2 g + 2 sinh g) 

l + ^^feF-nh^cosh 2 ^) ' 

( 2 sinh 6 f + 10 cosh 2 f sinh 4 f + | cosh 4 f sinh 2 ^ 
X = ^ 

1 4. 27A 2 
1 T 2m 2 M 2 „ 



gives 
(69) 

(70) 



27A 2 • u <*) h 2 *(*) 

2^Fm^P &mn — cosn — 
^ d^siBh*^ cosh (cosh 2 f 



By using ([3T]) and 



V 



^ cosn !a( cosh ^ +2 sinh^^ 



, we obtain 

27A 4 



\ (2 sinh 6 f + 10 cosh 2 f sinh 4 f + | cosh 4 f sinh 2 ^ 

<r(t) - 4M 2 ^ ^ „ 

1 + rajp smn — cosn — 



TV-72 sinh ^p- cosh 2 

eff 

lnh l im C o S h ^ ( cosh2 + 2 sinh ^) ) > 

1 4. - nVl *(*) h 2 *(«) 

\ 1 + 2m*M* ii l* S1Uh ~T CO&n — J 

27 A 2 . t(t) 3 t(t) 
■ sinh —7— cosh —7— , 

27A 4 



/ 27^ cinh t <!) 



z 2 



( (2 sinh* f + 10 cosh 2 sinh 4 ^ + § cosh 4 f 
-M 2 { ^ 

1 4. 27A 2 
1 ' 2m ! U a , 



a sinh 2_a N 

^^sinh^cosh 2 "^ 



(71) 
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27 A s KvnVi 2 ^W. cosh * ' 



cosh 2 l -Y- + 2 sinh 



2 7 A 2 



sinh co: 



2^ 



+4m 2 M 2 ff 



189A 2 .,*(*) l2 *(*) 
— — sinh —7^- cosh -^- L 
V I I 



(72) 



One may find £ as a function of t — £ by using the expression of a in (|39[) . Then in principle t is given as i = 
Substituting £ = t(£) into the expression U (t) in (|T2"|) . we can find the expression of C/ as a function of £, C/ = U (t (£)), 
which shows, by using the U in (|39[) . the expression of U((,). On the other hand, by comparing the expressions 
of V in (|39|) and ([70|) . we find V(ip). Then by following the procedure from (fTF)) to (fl"9|) . we get the expression 
of (R(-s\RV\e n (Vs -1 /)) and e„ f\/s _1 /))- Thus > thc ACDM universe can be realized 

without dark matter. This may suggest that the massive spin two particle might be a dark matter. In the same way, 
the reconstruction of F(R) bigravity realizing the given cosmological evolution may be done. 



V. SUMMARY 



In summary, we proposed a bigravity analogue of the F(R) gravity. Our formulation is based on recent ghost-free 
bigravity theory. The scalar fields are added in both metrics sectors of theory so that after corresponding conformal 
transformation the scalars become auxiliary ones. Integrating out auxiliary scalars, ghost-free F(R) bigravity follows. 
It turns out, however, that construction in terms of auxiliary scalars (i.e. when F(R) is given implicitly) is easier to 
work with. Cosmological equations of the theory under investigation are shown to be consistent. The cosmological 
reconstruction scheme is developed in detail. It is demonstrated that almost any evolution of physical universe may 
be realized while second metric solution which often could be flat space exists. The examples of cosmic acceleration 
which describe phantom, quintessence or ACDM universe are presented. The fact that ACDM universe may be 
realized without CDM indicates that massive graviton may play the role of dark matter. 

Of course, physical properties of F(R) theory under investigation as well as its other formulations should be further 
investigated. In this respect, note that it is difficult to get the explicit presentation of usual F(R) gravity which 
realizes arbitrary cosmological expansion since the reconstruction is made via the solution of the differential equation 
[20| . In case of F(R) bigravity, we can construct models directly in terms of the auxiliary scalar fields although it is 
more complicated to give an explicit form of F(R). 

We have not discussed the local tests of theory as well as the possibility to generate the fifth force which might 
not be neglected by experiments. We may construct a model which avoids such problems by using the Chameleon 
mechanism [2l| as in usual F(R) gravity (22j. An analysis by using the post- Newtonian parameter 7 was done in [23J. 
Such an analysis could be also applied to the models proposed in this paper. Moreover, the Vainshtcin mechanism 
Q might work to suppress the fifth force in general bigravity models. Furthermore, in case of the standard F(R) 
gravity it was proposed and studied Palatini formulation (Refs. [UH^ and references therein). Such formulation uses 
different variables set (connections) if compare with metric formulation. Formally, it may lead to the results which 
are not equivalent with the ones in metric approach. The investigation of massive bimetric F(R) gravity in terms of 
Palatini-like formulation looks an extremely interesting problem. For instance, does the ghost-free structure of theory 
survives in Palatini approach? This will be discussed elsewhere. 
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